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Starting in 2001, the MILC Collaboration began a large scale calculation of heavy-light meson decay constants 
using clover valence quarks on ensembles of three flavor configurations. For the coarse configurations, with a = 0.12 
fm, eight combinations of dynamical light and strange quarks have been analyzed. For the fine configurations, with 
a — 0.09 fm, three combinations of quark masses are studied. Since we last reported on this calculation, statistics 
have been increased on the fine ensembles, and, more importantly, a preliminary value for the perturbative 
renormalization of the axial- vector current has become available. Thus, results for fs , /b s , fn and fo s can, in 
principle, be calculated in MeV, in addition to decay-constant ratios that were calculated previously. 



1. INTRODUCTION 

Calculation of heavy-light meson decay con- 
stants is important for the extraction of a number 
of CKM matrix elements. For example, r(L>+ — » 
oc /!, U cs | 2 . As decay constants are fairly 
easy to determine in lattice calculations and they 
will be measured at BaBar, KEK and CLEO-c, 
they provide an excellent opportunity to verify 
the accuracy of our methods. 

We are extending a calculation of heavy-light 
meson decay constants with three flavors of dy- 
namical quarks that was begun in 2001 In 
this calculation, clover quarks are used for both 
the light and heavy valence quarks, the latter with 
the Fermilab interpretation 2 . A collaboration 
of Fermilab and MILC is now using light Asqtad 
quarks with the same dynamical gauge configura- 
tions to calculate decay constants. The newer cal- 
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culation |3] will allow better control of chiral ex- 
trapolations, as it is practical to reduce the light 
valence quark mass. However, for D s and B s we 
expect the two approaches to be complimentary. 

Since the last time we reported on this calcu- 
lation 0] , we have greatly increased the statistics 
on two gauge ensembles with a = 0.09 fm. (See 
Table 1.) Furthermore, preliminary results from 
a one-loop perturbative calculation of the axial- 
vector renormalization constant Za have been 
provided by El-Khadra, Nobes and Trottier 

Dynamical gauge configurations are generated 
using the Asqtad action [H]. Further details of 
confiuration generation may be found in Ref. 0. 
For each ensemble of dynamical quark configu- 
rations, we use five light and five heavy valence 
quark masses. The masses and decay constants 
are interpolated or extrapolated as explained be- 
low to get physically relevant values. The rela- 
tive scale is set through the heavy quark potential 
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0.031/0.031 
0.0124/0.031 
0.0062/0.031 


7.18 
7.11 
7.09 


496 
527 
592 


163 (163) 
242 (120) 
293 (48) 



Table 1 

Details of fine lattice configurations. In parenthe- 
ses are the numbers of configurations analyzed in 
Rcf. @]. 



[7] and the overall scale is set from bottomonium 
splittings. 

2. ANALYSIS OF RESULTS 

The analysis of the heavy-light decay constants 
involves a number of steps. On each ensemble 
studied, we: 

1. fit light pseudoscalar (PS) hadron propaga- 
tors to determine PS masses 

2. perform a quadratic chiral fit of squared PS 
masses to determine k c 

3. determine m s from the mass of ss pseu- 
doscalar state assuming a linear chiral mass 
relation 

4. fit heavy-light (HL) channels to determine 
their masses and decay amplitudes 

5. extrapolate or interpolate results in light 
quark mass to m U) d or m s , respectively (see 
Fig. 1). We see that the interpolation re- 
quired for the strange quark mass is well 
under control. However, in extrapolating 
to m u _d, we are expected to find chiral log- 
arithms 8 . With clover light quarks, it is a 
long extrapolation to the physical value of 
m Ut d and it is not possible to see the chi- 
ral logs. This difficulty is much improved 
with Asqtad light quarks (SJ ■ We could also 
consider an extrapolation of decay constant 
ratios, such as !b/U or (f B J 'f B )/(fK/ 7tt) 
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Figure 1. Chiral extrapolation (interpolation) for 
Sb (/b.) for /?= 7.09, kq — 0.119. 



6. after removal of perturbative logarithms, 
fit fQqy/MQq to a power series in 1/Mq 9 
and interpolate to B, B s , D and D s meson 
masses 

7. put the perturbative logarithm back and 
use the heavy-light axial-vector current 
renormalization constant to get the renor- 
malized decay-constant 

Just before the conference, we obtained the pre- 
liminary results of a perturbative calculation of 
the axial- vector renormalization constant Za |S] . 
As these results are preliminary, and our use of 
the results has not been as thoroughly checked as 
we would like, the following results are to be con- 
sidered preliminary. In particular, we have not 
yet tried tadpole improvement, which should help 
to determine the size of our systematic error. We 
note that when we reported on this calculation 
at Lattice 2002 4, , no perturbative or nonpertur- 
bative calculation of Za was available. We used 
an ad hoc procedure based on comparison of our 
improved action quenched results with the contin- 
uum limit of earlier calculations using the Wilson 
gauge action and Wilson or Clover quarks. This 
was explained in more detail in Ref. p^. 

After steps 1-7 are completed on each ensem- 
ble, we have a partially quenched result at a par- 
ticular value of dynamical /m p . We then plot 
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these results as a function of (m T /m p ) 2 to per- 
form a chiral extrapolation for the sea quarks. 
Before looking at the decay constant of a spe- 
cific meson, we note that if we plot the ratio of 
decay constants for two different mesons a good 
deal of the uncertainty from the renormalization 
constants and other systematic errors drops out. 
Figure 2 shows the ratio /b s //b of meson decay 
constants. In Fig. 3, we show fs s - 



3. FUTURE WORK 

We need to complete this analysis by includ- 
ing alternative cuts on the fits of meson propa- 
gators and alternative chiral extrapolations. We 
also must check the effect of tadpole improvement 
and see if we can employ a trick of Ref. JU| in 

which Zf = p^^Jz^Z 1 ^ where p% is computed 
perturbatively and both Zy's nonperturbatively. 
Although it would be possible to increase statis- 
tics on the fine configurations or include newer 
coarse ensembles (a = 0.12 fm), our more re- 
cent effort with Asqtad light quarks appears to 
be more promising. 
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